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Abstract 

We present a matrix model which is intended as a toy model of the gauge dual of an 
p) • AdS black hole. In particular, it captures the key property that at infinite correlators 

1^ , decay to zero on long time scales, while at finite N this cannot happen. The model 

consists of a harmonic oscillator in the adjoint which acts as a heat bath for a particle 
in the fundamental representation. The Schwinger-Dyson equation reduces to a closed 
^ . recursion relation, which we study by various analytical and numerical methods. We 

' discuss some implications for the information problem. 
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1 Introduction 



The black hole information paradox [T] is one of the great thought experiments in physics. 
Three decades of effort have made it clear that it has no trivial resolution, but indeed requires 
a modification of some central principle of physics [1] , although the focus has largely shifted 
from a breakdown of quantum purity to a holographic nonlocality of quantum gravity [2l[3]. 

The information paradox was largely responsible for the intense scrutiny given to the 
dynamical properties of black branes, which led to the discovery of gauge/gravity duality [1]. 
This duality in turn implies that information is not lost, because one can describe the 
formation and decay of a black hole within systems that have a well-defined dual description 
in an ordinary quantum framework^ 

This argument for information preservation is rather indirect. In order to calculate the 
black hole S-matrix, one must translate the initial infalling state into the dual field theory, 
evolve forward in the field theory variables, and translate back into the outgoing state of 
the Hawking radiation. It would be desirable to have a prescription entirely in terms of 
the bulk gravitational variables, since this is how the Hawking radiation, and the apparent 
information loss, is found. Indeed, the proposed answers to the question "Where does the 
argument for information loss break down?" are now as diverse as the answers to the original 
question "What happens to information thrown into a black hole?" We may hope that the 
attempt to answer this new question will be as fruitful as it was for the earlier one. 

It is notable that some features of the black hole persist even in the weakly coupled gauge 
theory, where the spacetime interpretation of the bulk breaks down. The continuation of the 
Hawking-Page transition [71 [8] to weak coupling has been studied extensively [9| [TOl fTTj [12]. 
It has also been argued that vestiges of the black hole singularity [I3] and of the information 
problem [H] survive at weak couplingEI 

Ref. [H] considers the information problem in the form presented in Ref. [15]. The low 
energy gravitational field theory description of AdS black holes shows quasinormal behavior, 
the exponential decay of correlations in time [191 120] • This description should be valid at 
large N and large 't Hooft parameter, where the curvature is small. In the dual field theory, 
exponential decay is possible at infinite A^, where the thermal field theory has an infinite 
number of states and can absorb an arbitrary amount of information. However, at large finite 

exponential decay can only persist until the correlations are of order e~'-^^^'^\ because the 
black hole has only a finite number of states and so the correlator is a sum of a finite number 

^In stating that the dual theory is well defined, we have in mind the superrenormalizable duals introduced 
in Ref. |5], as well as the BFSS Matrix Theory [Q, which is a quantum mechanical example of gauge/gravity 
duality. 

^Ref. [13], and our work, model the weakly coupled duals to the AdS^ x or Matrix Theory black holes. 
The weakly coupled limit of the BTZ black hole is also interesting, and the corresponding questions have 
been explored in Refs. [HI [H [HI [H] . 
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of exponentials. Reconciling this with theprediction of the low energy bulk field theory is a 
manifestation of the information paradoxjj 

Ref. [11] argues that signs of the quasinormal behavior can be seen at arbitrarily weak 
coupling, in that perturbation theory breaks down at long times no matter how small the 
coupling is. We would like to extend this, resumming the perturbation theory to see the 
explicit form of the late-time behavior. The graphs considered in Ref. [H] have a simple 
iterative structure that suggests such a resummation. However, we would like to find a situa- 
tion in which this resummation is systematic, in that it represents the full planar amplitude. 
This would provide a setting for discussing the breakdown of the large-iV approximation, 
which is dual to the loop expansion of the bulk theory. 

We identify a simple system with the desired property: a harmonic oscillator in the 
U{N) adjoint representation plus a harmonic oscillator in the U{N) fundamental, coupled 
through a trilinear interaction. In particular, it is sufficient to consider the limit of a single 
fundamental excitation in interaction with the adjoint oscillator. One can think of the adjoint 
as a heat bath coupled to the fundamental, and we study the decay of correlations of the 
fundamental field. 

This model may have a number of realizations; one can think of it as a reduced version of 
a DO-brane black hole with a DO probe, where the matrix variables are the 0-0 fields and the 
fundamentals are the 0-Oprobc fields [22], [23]. It might actually arise in some decoupling limit, 
although of course it would be a limit of large spacetime curvature. This simple model may 
also have applications outside of black hole physics, along the lines of the Caldeira-Leggett 
model [21]. 

In Sec. 2 we introduce the model. It has the same graphical structure as the 't Hooft 
model of two-dimensional QCD [25]. Unlike the 't Hooft model there are dynamical adjoints, 
but we achieve the simplified graphical structure of the 't Hooft model by stipulating that 
the adjoints have no self-interaction. We first consider zero temperature; we derive the 
Schwinger-Dyson equation, which reduces to a two-term recursion equation with respect to 
frequency. The singularities consist of poles on the real axis. We also analyze the quantum 
mechanics canonically, leading to a closed-form solution a zero temperature. 

In Sec. 3 we consider nonzero temperature. Again, the Schwinger-Dyson equation reduces 
to an iterative equation, but with three terms rather than two. In this case we can argue 
that the correlator for the fundamental field must approach zero at long times (in the planar 
limit). Numerically, we show that the decay is power law for small couplings and exponential 
for couplings that are sufficiently large; the power law regime is likely an artifact of our model 
that has no analog for the black hole. We develop briefiy the approximate solutions for small 

•^The relation between the information paradox and the discreteness of the spectrum is discussed further 
in Ref. 121]. 
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coupling and small mass. We also extend the model, by the introduction of decoupled sectors 
and the singlet constraint, so as to obtain a Hagedorn transition. 

In Sec. 4 we discuss some implications for the information problem. 

2 Zero temperature 
2.1 The model 

The fields are a Hermitian matrix Xij{t) and a complex vector (piit), with conjugate momenta 

[Xij,Uki] = iSiiSjk , [4>i,Tij] = iSij . (1) 

The Hamiltonian is 

H = ^Tr(n2) + ^Tr(X2) + n\l + gX/M)n + M''(t)\l + gX/M)<p . (2) 

In terms of the lowering operators for the fundamental and antifundamental, 

_ Tx\ - zM(f)i _ _ 71, - iM4 

this is (dropping a constant) 

H = -Tr(n2) + — Tr(X2) + M{a^a + a^a) + g{a^Xa + a^X^a) . (4) 

We do not impose the singlet constraint for now, and so there is no Hawking-Page transition. 
In Sec. 3.4 we will extend this model to include these features. 

There are several motivations that lead us to this model. The model of Ref. [H] is based 
on iteration of the basic graphical unit shown in Fig. 1. In our case, the basic process is 
just the upper half of this graph, above the dashed hne: a vector emitting and reabsorbing 
an adjoint. The doubling of the graph in Ref. [H] plays no essential role, so our model 
should have similar properties. Moreover, we will see that for us the iteration represents 
the full planar approximation, whereas in Ref. [H] it is just a partial summation, while the 
full planar summation would be much harder. For the purpose of systematic study of the 
breakdown of the 1/N approximation, it is useful that one can sum the full planar amplitude. 

In fact, a very similar model has already been studied in Refs. [221 [23] ^in approximation 
to the quantum mechanics of DO-branes. That model is more elaborate, in that the fields 
carry additional indices, and there are fermions, with supersymmetry. However, the basic 
graphical structure is the same; in particular the self-interaction of the adjoints is replaced 
with a quadratic potential, via a mean field approximation. In Refs. [221 123], the adjoints 
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Figure 1: The basic graphical unit studied in Ref. [14]. Iteration of this leads to breakdown 
of perturbation theory at long times. Our model iterates a basic unit which is just one side 
of this, above the dashed line. 

form a DO-brane black hole, and the fundamental is a string stretched from a probe DO- 
brane to the black hole. We will be studying the correlator in the one-fundamental sector. 
Essentially, we are looking at waves traveling on the stretched string, particularly as they 
fall into the black hole. (Ref. [26] considered a similar situation.) 

This system has no ground state, because the highest term is cubic in the fields. 
However, the Hamiltonian commutes with the number operators A*",^ = a^a and A^^ = a^a, 
and in each (A*",^, N^) sector there is a ground state. Defining 

H' = H + ciN^ + N^){N^ + N^-l) , (5) 

the eigenstates of H' are the same as those of H, and for sufficiently large c and M the 
ground state will be in the sector = = 0. We assume this henceforth. 

In fact, for our purposes we need study only the behavior of a single particle in the 
fundamental represention, in interaction with the matrix heat bath. To isolate this we will 
take the splitting M for the fundamental oscillator to be large compared to all other scales, 
in particular the temperature, and we focus on the observable 

^iMit-t') (ja,{t)a]{t'))^ = 5,,G{T,t- t') . (6) 

Note that t and t' are Lorentzian, and that in the present section we are interested in 
temperature T = 0. By construction, the stabilizing term ([5]) vanishes in the relevant 
sectors A^,^ = 0, 1, A''^ = 0, and so we can calculate with the original H (j4]). Including the 
phase factor in the correlator, the dependence on M drops out in the large-M limit. Thus 
there are essentially two parameters at zero temperature, m, and g with units of m^/^. In 
terms of the analogous brane system, taking M to be large means that the probe brane is 
far from the black hole. 
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Figure 2: Schwinger-Dyson equation for planar contributions to G{lj) (propagator with 
shaded rectangle) in terms of Go{u!) and Ko{uj). 

2.2 The Schwinger-Dyson equation 

At zero temperature, the sum of all planar contributions to the correlator (j6]) is given by the 
Schwinger-Dyson equation 

—G{J)IU{uj - uj') , (7) 



where X = g'^N and 



G.iu) = , Kn(u) 



iu + ie ' " ' ' Lu"^ — rn? + ie 



This is shown graphically in Fig. 2. This has the same form as in 2-D QCD [25] because 
the index structure of the interaction is the same. The integral can be carried out. Because 
a annihilates the vacuum, G{t) vanishes for t < and so G{uj) is nonsingular in the upper 
half-plane. Also, because the coupling g has mass dimension 3/2 we can assume that G 
has its free behavior i/uj aX high frequency. We can then close the contour in the upper 
half-plane and evaluate the residue ai uj' = uo — m + ie io obtain 

G{uj)='-[l~^G{uj)G{uj-m)^ . (9) 

We omit the ie from this equation, with the understanding that the correlator is to be 
evaluated infinitesimally above the real axis. 

In Sec. 2.3 we will solve this in closed form, but first let us study it using a variety of 
analytic and numerical approaches. To get our bearings, let us consider first the limit m — 
with 2A/m = z/^ fixed. Eq. ([9]) becomes the algebraic equation 

z/^G^M -4zcjG'(cu) -4 = , (10) 
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and so 

2i 



(11) 



UJ + \/ iJ^ — V 



On the physical sheet, the square root approaches lo at long distance. The u; = pole has 
been broadened into a branch cut of width Iv. 

This has a simple interpretation. We are taking the harmonic oscillator frequency m 
to 0, so the matrix X is essentially a static variable with an eigenvalue distribution given 
by the Wigner semi-circle law. The propagator (ITT!) is just the free propagator with mass 
/i = (yfX, averaged over a semi-circle eigenvalue distribution for X of width ^j2N/m. The 
branch cuts in uj translate into the asymptotic power law decay e^*''*^"^''^. However, this 
power law decay is unrelated to the quasinormal behavior that we seek: it originates from 
the noncompactness of X due to the vanishing of the potential at m = 0, rather than the 
large- limit. 

Writing Eq. ([9]) as a recursion relation 



G{uj) 4 



G{LO-m)~iuj. (12) 



gives an efficient way to calculate numerically, beginning with the asymptotic behavior G{uj) ~ 
i/{uj + ie). In order this procedure to work, we must require that the recursion be stable at 
large — otherwise, subasymptotic terms could grow to become significant. Suppose that 
there is a solution G^{uj), and we consider a perturbation G^iu) + 7(0;). Then 

7M = -^G^M7(^-m) . (13) 

The recursion is stable towards increasing oj when \u'^G1{ijj)\ < 4, and stable towards de- 
creasing u when the inequality is reversed. Since G{uj) = 0{uj~^) at large |a;|, the recursion 
is stable there in the direction of increasing u. Given perfect numerical precision, we could 
start at very large negative u and use the convergence to bring us very close to a solution. 
Even if the recursion became unstable for some intermediate range, we could make the initial 
error as small as desired. In practice, a highly unstable recursion would lead to a numerical 
errors. This is an potentially an issue at small m, where the recursion requires many steps. 
In fact, using the m = solution (fTTl) one finds that the recursion to the right is no worse 
than neutrally stable in this limit. 

Before performing the numerics, we can anticipate that the branch cut found at m = 
must break up into poles. Since there is no branch cut at large u, none can appear through 
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the recursion. However, whenever the right-hand side of Eq. (1121) vanishes, there will be a 
pole in G{uj). Generically this will happen at isolated points on the real axis, because the 
right-hand side is purely imaginary there. The numerical integration verifies this picture: 
the branch cut breaks up into poles as m is turned on, and the poles move further apart as 
m is increased with i/ fixed. Zero temperature results for the real part of G{uj) are shown in 
Fig. 3 for u = 1, m = 0.05 and in Fig. 4a for u = 1, m = 0.80. 

There is a temptation at small m to approximate the recursion ( |T2l) by a differential 
equation, but this does not seem to be useful: the differential equation is less stable than 
the recursion relation. They agree on smooth configurations, but have different e^^'^^"^^ 
instabilities. 

2.3 Canonical calculation 

The intermediate states in Fig. 2 consist of a single (p excitation plus any number of X 
excitations, with the indices contracted: 

\j,r) = N-^/'aliA^%\v) , r>0, (14) 

where \v) is the free ground state and A]- = (Iljj -|- iraXij) / \/2m is the raising operator. The 



states |j, r) are orthonormal at large A^. Then 

. r-l 



{H-M)\j,r) = mr\j,r) + -\j,r + l)-'-\j,r-l)--Y,N^'~''\A^%k\j,r-l-l) . (15) 

1=1 

The last term has norm of order A^^^ and so drops out in the planar limit. The noncompact- 
ness in r suggests that this model may show quasinormal behavior at finite temperature. 
Defining 

oo 

|j,^) = J]^,|j,r) , (16) 

we have the eigenvalue condition 

{uj - mr)^pr = ^{A-1 - A+i) , ^^-1 = 0. (17) 
This is essentially the Bessel recursion relation. Specifically, 

V'r = i"" Jr-ulm{v/m) . (18) 

Here Nr-uj/miy /f^) would satisfy the same recursion relation, but only the solution (fTSl) 
is normalizable. The eigenvalue condition J~i^uj/m{i^ /^) = determines the poles in the 
correlator. With this clue, we can find the closed-form solution to the Schwinger-Dyson 
equation ([2]), 

n( ^ 2z J_^/^(z//m) 
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ReG(a>f0.10 



b) 



-1 1 CO 

Figure 3: a) The real part of G{lj) for i' = 1, m = 0.05, evaluated 0.01 units above the real 
axis to give the delta functions finite width, b) The same function evaluated 0.1 units above 
the real axis: the poles merge into an approximate semicircle distribution. 
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3 Finite temperature 



3.1 The Schwinger-Dyson equation 

Now consider the system at nonzero temperature. To study real-time thermal correlators 
one generally needs the doubled integration contour of the Schwinger-Keldysh formalism [271 
EHlEn]. However, our situation simplifies. First, we are assuming that M is large compared 
to the temperature, so there are no <^ excitations in thermal equilibrium. Second, the X 
fields have no self-interaction, so the thermal ensemble is freefl In this case the real-time 
formalism reduces to replacing the X propagator with the free thermal propagator, 

i [ 1 e""*/^ \ 
^o(r, u;) = zwT -2 ^ • (20) 

The Schwinger-Dyson equation is changed only by the use of this propagator, 

—G{T, uj')Ko(T, uj-uj') . (21) 

-oo 

The time-ordered correlator G(T, t) still vanishes at t < (because M ^ T) and so we can 
again close the contour in the upper half-plane and pick up only the poles of Ko{T, uj). This 
gives 

G(T, ..) = ^ 1 1 - ^^^-^^-^G(T, uj) [g(T, u;-m) + e-"^/^G{T, ^ + m)] | , (22) 

or 

G(r, u-m)- -j^-— + e-/^G(T, u; + m) = ^ (23) 

with 4 = 1^7(1-6-"^/^). 

Although this is similar to the zero-temperature Schwinger-Dyson equation, the behavior 
of its solutions is very different. At zero temperature the singularities of G{uj) are poles on 
the real axis. At nonzero temperature such poles are impossible. 

To see this, note first the spectral representation, 

G(T,t-t') =ArTr(e-^/V(t)aI(t')) = AT^ |(A|ai|5)p0(t - t')e-''^/'^-^(*-*')(^^-^-^) , 

A,B 

(24) 

where Af^^ = Tre"^''^. The Fourier transform is 

GiT,.)^^^y:^-:^^^Mhm 



to — n 



(25) 



■^Note that in any event the backreaction of the fundementals on the adjoints is suppressed in TV. 
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We have introduced the nonnegative spectral density 

F{ii)=UY,e-''^'^ma,\B)\^5{ii-EB + EA) . (26) 

A,B 

Now, suppose that there is a pole in G(T, uj) at some ujq. Then the first term on the LHS of 
Eq. (l23l) has a pole aX uj = ujq + m. This must be cancelled either by a pole in the last term 
with a residue of the opposite sign, or by a pole in the second term from a zero of G(T, uj) at 
UJ = ujq + m. The spectral representation forbids a negative residue, so the zero must exist. 
Similarly we can conclude that there is a zero at ujq — m. It then follows that all three terms 
on the LHS vanish at uj = ujq, which is a contradiction except possibly for a pole at cuq = 
(which is indeed present in the free theory). 

The absence of poles, and so of delta-functions in F{fi), immediately implies that the 
correlator goes to zero asymptotically in time: this model has the planar behavior that we 
seek. If the zero temperature poles separate into branch cuts on the real axis there will 
be power law falloff; if the singularities drop below the real axis onto the second sheet the 
decay will be exponential. We find numerically that at small nonzero temperature the poles 
widen into branch cuts, while at higher temperature the cuts merge and the singularities 
drop below the axis. We will present these results in Sec. 3.2, after some further analytic 
discussion. 

The recursion relation again becomes algebraic in the limit m — ^ 0, now with z/|. fixed, 
(1 + e-"'/^)u^G\T, u) - AiuG{T, u) - 4 = . (27) 

The solution is 

2i 



G(T, uj) = ^^^^.m/ry^ I - " \/-^ " (1 + e-/^)4 ) ■ (28) 



The logic is the same as at zero temperature, with the eigenvalue distribution thermally 
broadened. 

The three term recursion relation fl23|) is unstable in both directions, so its solution is 
less constrained than at zero temperature. We can derive some useful results from the real 
part ReG(T,uj) = 7iF{uj). Then 

F{uJ - m) - F{cj) + e-'^'^Fiuj + m) = . (29) 

First, if the spectral density has support in any segment of the real axis then it has support 
in every segment translated by a multiple of m. Thus, when the poles spread into branch 
cuts, each cut is accompanied by an unbounded series of additional cuts. This does not 
contradict the known asymptotic behavior, because the magnitude of F{uj) goes to zero at 
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large uj. The form of the falloff follows from the fact that coefficient of the middle term in 
Eq. ( !29l) becomes large asympotically, 4co'^/ v\. The recursion relation is then dominated by 
two terms, 



F{uj -m)lF{uj) 
F{u + m)/F{uj) 



= e 




— CO 



(30) 



so the spectral density behaves asymptotically as |u;| '^d'^D. 

3.2 Numerical results 

The instability of the recursion relation makes numerical calculation challenging. One ap- 
proach would be to fine tune the initial condition. However, we have had more success by a 
different approach, solving the stable zero temperature recursion relation and then solving 
the coupled differential equations obtained by differentiating the recursion relation (l23l) with 
respect to T at fixed zz-p. Figs. 4 show F{uj) for a series of temperatures from zero to infinity, 
with ut = I and m = 0.8. 

At zero temperature (a) this is a set of delta functions, meaning poles in G{u!). At low 
temperature (b) the poles become short branch cuts, and additional cuts appear, shifted 
by multiples of m. As the temperature is increased (c) the cuts lengthen and begin to 
merge, but there are still gaps. At higher temperatures (d) the cuts have completely merged 
and F{uj) is everywhere positive and smooth: the singularities have moved onto the second 
sheet. The recursion relation and its solution have a sensible infinite temperature limit (e). 
In the infinite temperature limit the equation is symmetric under G{uj) G{—uj*)*, and 
this symmetry provides a check on the numerical approach. We have verified numerically 
that the high-temperature solution is smooth out to |a;| = 4, and the asymptotics (130|) imply 
that nonanalyticity cannot arise at large u if not present at smaller values. 

The Fourier transform of Fig. 4e is shown in Fig. 5. Two exponentials are evident. The 
amplitude and decay rate of the first few points matches the central peak in Fig. 4e. The 
long-lived exponential arises from the closest pole to the real axis, which has a residue smaller 
by a factor of order 10"^. This appears to be associated with the last branch cuts to merge, 
and with the kink near u = 0.3 in Fig. 4d (examination with greater resolution verifies that 
the function in Fig. 4d is smooth at this point). 

3.3 Some approximations 

Thus, our simple model has the behavior that we seek. It would be good to have an approxi- 
mate analytic treatment of the quasinormal behavior. This is difficult because the model has 
three energy scales, m, ut, and T, and the quasinormal behavior disappears if any are set 
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Figure 4: The real part of G{uj) for vj- = 1, 



m 



0.80, and various values of y 



-m/T 



The vertical axis is rescaled at each temperature for best visibility (the actual area under 
the curve is vr at all temperatures), and at zero temperature ut is taken slightly above the 
real axis for the same reason. 
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Figure 5: The logarithm of the real time infinite temperature correlator, In |G(t)|, for vt = 1, 
m = 0.8, and T = oo. For clarity, time differences of 27v/m are displayed: on shorter intervals 
the correlator shows strong oscillations due to interference of singularities spaced by multiples 
of m. 

to zero|§ Note that in Fig. 4 we have taken z/^ = 1, m = 0.8, near the middle of parameter 
space, because the numerics are cleanest there. 

We can develop an approximate solution when any of the parameters is small; here we 
note a few features of two such approximations. 



3.3.1 Small A (small i^t) 

The ordinary perturbation theory in A is singular here at long times [2]. In our model, this 
shows up as the fact that each additional order of perturbation theory brings additional poles 
in G{T,uj), even within the IPl part. However, we can obtain an improved approximation 
from what we know of the solution. We expect that the pole at u; = will turn into a short 
branch cut, and that additional branch cuts will develop translated by multiples of m. From 
the structure of the recursion relation, the spectral weight in the branch cut at u = mk will 
is of order gf"^'*'!. Again, z/|, = 2A/m(l -e-™/^). To see the branch cut we focus on <^ m. 



■''There is power law decay at m = 0, but as we have noted this is not associated with large N. 
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for which the recursion relation imphes 



2 ' 
T 



z/; 



'TG{T,uj + m) 
^ --I^G{T,uj) = -^ . (31) 



iyiG{T,u-m) K 

In the second and third hnes we have have used <^ m on the RHS, and have dropped one 
term on the LHS because G{T, uj±2m) ^ G(T, uj). Using these we can ehminate G{T, uj±m) 
to obtain the cubic equation 

iujyG{T, uj)'^ + (-1 + y ~ - CUT + yuT)G{T, tof + ir(2 - 2t/ + ut)G{T, uj)+t^ = , (32) 

where y = e~"^^'^ and r = Am/u^. For simphcity consider the hmit y ^ 1, where 

iuG{T, - G{T, uy + iLJT^G{T, u) + = . (33) 

At u; = 0, G{T, cj) = r is purely real. However, the discriminant for this equation vanishes 
at u;^ = (11 + 5^/^)/2r^, beyond which the solution is purely imaginary. This determines the 
branch cut width 0(r~^) = 0(A/m^[l —y])- 

These results have been confirmed numerically. Note that the spectral density implied 
by the cubic equation is more complicated than the semicircle law that might have been 
expected for a short cut|§ The cuts at ± m are also determined by the relations (ISTIl and 
are smaller by 0(A). The further cuts are given by the two-term recursion fl30l) and are down 
by A''^' as expected. One can systematically carry this small-A approximation to higher order 
and no singularities arise, so we have found the quantitative small-A behavior. The branch 
points imply power law decay, so there is again a contradiction with the finite- behavior. 



3.3.2 Small m 

At small m the recursion relation goes over naively to the differential equation 

m(l - e-"^/^)G'(T, u)G{T, u)' = (1 + e-'"/^)G(T, - ^G(T, oo) - \ , (34) 

where the prime denotes an cu-derivative. We are holding m/T fixed. If T were held fixed 
the LHS would be of order and we would need to keep also a second derivative term, 
but the conclusions would be similar. The derivative term is a perturbation, but a singular 

^The discriminant also vanishes at = (11 — 5'^/^)/2t^, so there is another branch cut not far below the 
real axis. 
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one [30]: it can become large in regions of large gradient. Again it is useful to study the 
stability near a solution G^,, expanding G = + 7: 

^(1 _ e— /^)a7' = - e-'"/^)G:7 + 2(1 + e-'"/^)a7 - ^7 • (35) 

Assuming that the solution G* is close to the algebraic solution (1281) . one finds that the 
solution is stable toward increasing u for u; < —u and for > where ±z/ = iz/yVl + e~"^/^ 
are the branch points of the m = solution, and stable toward decreasing uj ioi —0 < uj < v . 

Thus, starting from the known asymptotic behavior at large negative u;, the differential 
equation has a unique solution which is close to the m = solution for uj < —v and then 
blows up rapidly, as e^^ . Similarly by starting at = i> and integrating in both directions 
one obtains a solution that is close to the m = solution for uj > —v and blows up beyond 
this point. These solutions clearly do not match onto one another as a single solution to 
the differential equation. The point is that near the branch point lo = — z> the gradients 
become of order 1/m and we must use the original discrete form of the equation. The full 
solution requires matching this 'boundary layer' solution with the smooth outer solutions 
that we have found. The details are an interesting direction for future work. In particular, 
the singularities will be closest to the real axis in the boundary layer, where the variation is 
the most rapid, and so this part of the solution dominates the long-time behavior. 



3.4 The Hawking-Page transition 

In the model thus far, we have not imposed the singlet constraint that generally is present in 
quantum mechanical realizations of gauge/gravity duality. In consequence, there is no phase 
transition as the temperature is varied. These features can be restored without altering the 
earlier results, by introducing some decoupled sectors. First, to form a singlet we need an 
antifundamental excitation, which we can obtain from an additional decoupled oscillator. It 
would be more natural to excite an antifundamental excitation of the original (p oscillator, but 
then we would also need to include ladder graphs as in Ref. [25]; this may be an interesting 
direction for future work. 

The thermal phase transition with large- oscillators has been studied in Ref. [31] and in 
more detail in Ref. [12] . We will review this very briefly. The singlet constraint is enforced by 
integrating over a Wilson- Polyakov line U in the Euclidean time direction. At fixed diagonal 
U the adjoint thermal propagator is 

/oo 
dte''''{TX,jit)Xkim=5uS,k 
-00 

00 -X 
+^ E [ie-"'^''U,,Ur^'m^ + H + {e-/^U,,U^')-6{u - m)] . (36) 



{ i 
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The sum represents paths that wind around the Euchdean time direction ±n times, each 
winding picking up a phase from the Wilson- Polyakov hne. All adjoint propagators in a 
graph are evaluated with the same U, which is then integrated over. The propagator is 
unaffected because M is large: paths that wind around the Euclidean time direction do not 
contribute. 

We extend the previous model by the addition of an adjoint oscillator of frequency m'. 
This is decoupled from the other fields except through the singlet constraint, that is, through 
its coupling to U. The phase transition occurs when [12] 

As m' — ^> 0, Tc — > and as m' — > oo, Tc ^ oo, so we can adjust the transition to take place 
at any temperature. 

At T < Tc, the eigenvalues Ua are uniformly distributed on the unit circle, so that 
Xlt ^ii ~ finite n. It then follows that the fZ-dependent terms from Eq. (136!) drop 

out in every trace at leading order in A^, so the correlator reduces to its zero temperature 
form (see for example Refs. [321 [33]) Thus there is no quasinormal behavior. 

At T > Tc the eigenvalues become nonuniform, and at T ^ Tc they are concentrated 
around the identity. For T ^ m' we can replace U with the identity, and the propagator fl36l) 
reduces to its form ( !20l) without the singlet constraint. We then recover the quasinormal 
behavior found earlier. 

4 The information paradox 

Now let us return to the information paradox, and discuss what we might learn. Our model 
does not capture all aspects of the paradox, spacetime locality in particular: like the weakly 
coupled gauge theory, there is no large bulk dual and no notion of spacetime locality. At 
best, the model represents a bulk theory in a large curvature limit. 

What the model does preserve is the large- A^ structure, the existence of dissipative behav- 
ior in the planar limit that does not survive at finite A^. The 1/A^ expansion in gauge/gravity 
duality is dual to the loop expansion in quantum gravity. It is a crucial question, how the 
preservation of information is manifested in this expansion. Whatever form the answer takes, 
there should be some parallel in our model. 

For example, there have been many proposals over the years that evidence for informa- 
tion restoration can be seen in the breakdown of the gravitational loop expansion even at 
low orders; see Ref. [31] for a recent discussion along these lines. A rather different proposal 
is that the restoration arises from nonperturbative effects in the gravitational loop expan- 
sion, additional saddle points that were omitted in the original argument [151 IM] . These 
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two approaches would correspond respectively to an order-by-order breakdown of the 
expansion, and to the contribution of e"*^^^^-* effects. 

The saddle point proposal has been analyzed critically in Refs. [HI [171 [36l [371 [33 [39] . We 
can paraphrase their key argument as follows. In the planar approximation, in the regime 
where the singularities lie below the real axis, the spectral density F{fi) = ReG{fi)/7i has 
support on the whole real axis. At finite A^, this must break up into poles with a typical 
spacing e~'^^^^^ If we consider an observable corresponding to the convolution of G{fi) with 
some smooth function, the effect is only of order e~'^^^^\ and so looks as though it might be 
captured by a saddle point contribution [151 [35] • However, if we measure G{fi) at a precise 
value of /i, the effect is of order one (or larger), which cannot be captured by such a saddle 
point. 

A Euclidean saddle point gives a contribution to the density matrix, so with two such 
saddles we have 

P = 1 + e^i-g^ ' ^^^^ 
with pi and p2 normalized density matrices and 5*2 — 5*1 of order G^^ ~ A^^. This implies 
for the correlator 

p = . (39) 

Observing the order one effect requires measurements on a time scale of order e*^'-^^-', which 
suggests the possibility of a compensating A^-dependence|ll However, the extra saddle points 
do not have the necessary enhancement: we need that the total residue of the poles in -F(p) 
be 1. The dominant saddle gives a continuum with total weig ht 1 - e-°(^ ), while the 
secondary saddles may give poles but with total weight e~'-^^^^\ 

Our discussion suggests a further problem with this idea. Namely, the difference between 
the planar and finite- A^ behavior already shows up in the basic model presented in Sec. 2.1, 
without the singlet constraint and without the Wilson line variable U. The additional saddle 
discussed in Refs. [IHl [35] appears in the [/-integration. 

We should be alert to possible artifacts of our model. In particular, the fact that we find 
power law decay at small g for arbitrarily large temperatures is likely connected with the 
fact that the fields in our adjoint heat bath are free, so states do not mix as completely as 
they should. We expect that in the weakly coupled gauge theory (at temperatures above the 
Hawking-Page transition) there will be exponential decay. However, the power law decay is 
already in contradiction with the finite- A^ behavior, so for our purposes we can regard this 
as quasinormal behavior. 



^For black holes that decay, the time scale is very much shorter. However, in this case one must observe 
of order particles, so again there is the possibility of an offsetting factor. 
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Because we do not have a bulk spacetime, we cannot separate stringy physics from 
gravitational field theory. That is, we capture the gravitational loop expansion but not the 
a' expansion. Thus we cannot directly investigate proposals such as those in Refs. [101 1^ . 
which relate information recovery to specifically stringy physics. 

Our discussion has been entirely in the field theory language. The dual language would 
correspond to working with U{N) invariants such as products of A and , either traced 
or in bilinears with a and a^; it would be interesting to develop this further. In this form, 
finiteness of shows up as relations between these invariants; for example, Tr(y4''") can be 
expressed in terms of lower traces for k > N. One does not expect these relations to be 
visible in perturbation theory in since for any given k they turn on abruptly at a finite 
value of Rather, they seem to imply that the closed string bulk variable are simply not 
good variables for nonperturbative gravity. 

This reduction of the number of independent variables is the 'stringy exclusion principle' 
[42] . It has been related to the growth of the size of objects at high energy [43]. The 
information paradox also implies another reduction in the size of the Hilbert space, the 
'black hole complementarity principle' [H]. For example, in an eternal AdS black hole 
there is an infinite number of infalling modes at the horizon; however, only a Bekenstein- 
Hawking entropy's worth can be independent. The idea that black hole complementarity is 
a consequence of the stringy exclusion principle is implicit in various places; the formulation 
of the information paradox in Ref. [15] makes it particularly clear. 

Growth of string states near the horizon was discussed in Refs. |l5l HB] . Unlike the case 
studied in Ref. [13], where the string blows up into a spherical D3-brane, at the horizon it 
grows into a randomly walking long string. In the gauge theory, the former corresponds to 
the trace of a large power of a single field, while the latter suggests the more generic case of 
the trace of a random sequence of fields. 

D-branes are related to the (2nioop)! growth of perturbation theory [47], versus nioop! 
in field theory, and to the associated e"^/^" effects (though they cannot be the only such 
effects [IH])- Given the connection between the stringy exclusion principle and D-branes, 
we are led to conjecture that the breakdown of low energy effective field theory in the 
neighborhood of a black hole is manifested by a (2r2ioop)! growth of perturbation theory, and 
by e"'^'-^/^-' nonperturbative effects, as compared to e~'^^^^'^^ contributions from field theory 
saddles. This conjecture is just based on analogy; we do not have any specific scenario for 
how these large high order amplitudes appear, and for how the e"*^*^^/^^ effects become 
order one in certain observables. 

Lastly, we note an interesting recent paper [19] on the time scale for black hole information 
return. The arguments of this paper rest on certain assumptions about the thermalization 
process. It may be possible to investigate these in our model, or some extension of it. 
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